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Introduction 

. 0.1 Adiabatic limit in quantum mechanics 

> I 

^ ! The problem of adiabatic limit in quantum mechanics with avoided eigenvalues 

crossings will be the basic example in our paper. 

Let us describe it in a precise way. We start with a smooth family of N x N 

^ Hermitian matrices A{t) with a <t < b and consider the following linear system 

. of differential equations: 

O ■ hdX ^ , ^ 

^ \ Solving this equation gives an unitary map Sh defined by Sh{X{a)) = X{b). We 

O ■ call it the scattering matrix of the problem. Adiabatic Theorems describe the 

asymptotic behaviour of Sh when /i — )■ 0. 



1. 



, , _. The clasical adiabatic theorem concerns the case where the eigenval- 

^ ■ ues of A{t) satisfy, for all t: 



Ai(t) < A2(t) < ■■■ < Aiv(t) . (2) 

In this case, if we start with Xj{a) an eigenvector of A{a) with eigen- 
value Aj(a), we have Xj{b) = exp{i J^^ Xj{s)ds)Yj{b) + 0{h) where Yj{b) 
is obtained by parallel transporting Xj{a) along [a,b] in the eigenbundle 
Lj{t) = keT{A(t) — Xj{t)) w.r. to the geometric (or Berry) connection de- 
fined by 

^d/dtY{t) = proji^(t) — 
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2. Another case which is well known [H |3l [131 ttS]) is the case where the 
eigenvalues of A{t) cross transversally. Then the previous results remains 
true with a less good remainder term 0{y/h) but we need to label the 
eigenbundle by following smoothly the eigenvalues at the crossings points. 
In other words, there is some exchange between Lj and Lj±i at each crossing 
point. 

So the scattering matrix remains diagonal but with some relabeling of the 
indices j = 1, ■ ■ ■ , N. 

We want to consider the bifurcation between the 2 situations. For that pur- 
pose, we allow A{t) to depend smoothly of a (i— dimensional parameter /i close 
to in R'^. We will assume that Ao(t) admits eigenvalues crossings transversally. 
For /i small, will have avoided eigenvalue crossings. We will derive in this pa- 
per a way to get the asymptotic behaviour of S^,h w.r. to both small parameter 
/i and h. 

Basically the result uses 2 parts : 

1. The local situation where a suitable Landau-Zener formula can be used 
describing the local 2x2 scattering matrix (see [T5| |8]). 

2. The global problem where we have to take into account interferences given 
by what we called Bohr-Sommerfeld phases. In the case of several cross- 
ings (analytic case), using Stokes lines, they are many works where the 
interferences effects are exponentially small [TB| [TS] . 

0.2 Born-Oppenheimer approximation 

Another important example comes from the Born-Oppenheimer approximation. 
We consider, on the real line a Schrdinger equation with matrix valued potential 
as follows: 

H = ~h^^^+V{x) (3) 

where V is an Hermitian matrix depending smoothly on x G M. We are interested 
in the eigenvalue equation {H — E)u{x) = with -u : M — t- C". If there is no 
degenerate eigenvalues of V{x) smaller than E, then eigenf unctions can be foud 
using a WKB-Maslov Ansatz in a very similar way to the scalar case. It is the 
generic situation for a single V [HI EH])- The most precise results are shown in 
|12j . However, if V depends on some external parameters /i close to in R'^, it is 
possible that, in a stable way, V^(a;) admits degenerate eigenvalues. The problem 
we adress in this case is the description of Bohr-Sommerfeld rules uniformly w.r. 
to /i. 
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0.3 The general setting 

Both examples can be put together as follows: we consider a self-adjoint semi- 
classical system of (pseudo-) differential equations of order with N complex 
valued unknown functions U = (f/i, ■ ■ ■ , Ujy), 

HU = 0(/i°°) (4) 

on the real line. Here H = (-^j,j)i<j,j<Ar is a matrix of semi-classical (pseudo- 
)differential operators of order in 1 variable x with {HijY = Hij. Viewing the 
revious ewmples, it will be important to consider the case where H = depends 
smoothly on a germ of parameters /i G (M"^, 0). In our previous papers [5l [6], we 
derived normal forms near the eigenvalues crossings which allow to compute a 
local scattering matrix including the Landau-Zener amplitude. The goal of this 
paper is to compute global objects including interferences effects. The general 
picture is already provided by the study of the scalar case j9j from which we 
know that we need to define ad hoc Bohr-Sommerfeld phases. 

The general terminology is the same as in [5] and [6] , but in the present paper, 
our phase space will always be 2 dimensional: 

ffl,, : n ^ Herm(C^) , 

the (matrix valued) principal symbol, is the dispersion matrix, and = p^^(O) 
with = det{H^^^J the dispersion relation. 

We first recall the local normal form as derived in our previous papers |5l |6] 
and we solve it. After that, we come to the new part which consists in deriving 
global objects in the spirit of [9]. 

We will need another piece of information which we call Bohr-Sommerfeld 
phase; let us take any simple cycle c (with singular vertices Zi, - ■ ■ ,Zj,--- ,Zp) of 
the dispersion relation Cq. We will associate to c a real valued symbol S/i(/i) ~ 
J2^o Sjif^)^'' were the S/s are formal power series in /x. 5*0 is a purely classical 
object which involves regularized action integrals. 5*1(0) is computed using the 
transport equation which is smooth along the edges [zj, zj+i] (Berry phases) and 
singular Maslov indices. From the Bohr-Sommerfeld phases we recover the global 
objects mod 0{h°°). 

1 The local normal form 

Let us recall the following result from [HI [5l [6] (see also [8]): 

Theorem 1 Let us assume that the function po{x,C,) = det(ifciass,/x=o) admits at 
the point zq G T*M a non degenerated critical point of Morse index 1 (also called 
hyperbolic critical point, because the Hamiltonian vector field of pq is hyperbolic 
at the singular point zq) and with critical value po{zo) = 0. 
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Then, we can find the following objects which depends smoothly on jj, close 
enough to 0; 

• A smooth family of germs of canonical transformations Xfi '■ {T*M., 0) — > 
{T*R, zo) such that 

with Cfj, an invertible germ function and 70 a germ of > function of fi 
satisfying 7o(0) = 0. Moreover, the Taylor expansion of '-/q is unique. 

• A smooth family of unitary FIO 's associated to Xfi o-nd N x N matrix 
of^DO 'sA^ 

so that, we have the following normal form ( called the Landau-Zener normal iorm) 
near Zq: 

^*ij.U^Hi^iUi^iA^ = 

with D = j-^, a{fi, h) T,^Qaj{fi)h^ a symbol and Q is elliptic. 
Moreover, Qq is a complex valued function of /i which satisfies 

|aor(/i) = 7o(/i) , 

and we have 

Remark 1 In JM/, Theorem U\ is proved under the following transversality hy- 
pothesis: 

{^) ii W C Herm(C^) is the submanifold defined by dimkeriJ = 2, we 
assumed there that {fi,z) — )■ i/ciass,^t(^) is transversal to W at the point 
(0,^0). 

This hypothesis can be restored using more parameters, so that Theorem U\ is 
also correct. For simplicity, 

we will assume that hypothesis (*) holds true in what follows. 

It implies that /i — )■ ao(yu) is a submersion from (M'^, 0) onto (C,0). We will 
denote by Z = ^(0). Z is a smooth germ of codimension 2 manifold of (M.^, 0). 
Formal expansions w.r. to the parameter fi mean formal expansions along Z . 

Remark 2 Contrary to the scalar case, there is no arbitrary choice concerning 
the images of the half axes = 0, x > 0}, ... by x- The smooth arcs of the 
dispersion relation are oriented in the following way: there is a change of the 
Morse index of the quadratic form associated to i/ciass from m to m ± 1 while 
crossing these arcs. The sign of this change is preserved by the gauge transform 
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Figure 1: the jumps of the Morse index of the dispersion matrix 



which acts directly on the previous quadratic form by an invertible linear change 
of variable. We choose to orient the arcs so that the Morse index is bigger on the 
right than on the left of the path. 

Remark 3 The symbol a is not uniquely defined because a diagonal unitary 
gauge transform preserves the normal form while changing a by some phase shift 
exp{iip{h)) . Its modulus 7(/i, h) = \a{fi, h)\'^ is uniquely defined from the Landau- 
Zener coefficient given in Equation 

The matrix A is defined up to matrices which will change the phase of a. 
More precisely, if Aq is the principal symbol of A at the crossing point, the 
only prescription is that Aq is a map from to which sends © into 
keTHo{zo) and satisfies Aq ((C © 0) © 0) = Ei and Aq ((0 © C) © 0) = ^2 where 
Ej = hm^^^g^^gAA^^okeriJo with Ai = xo({^ = 0}) and A2 = Xoi{x = 0}). The 
choice o/(Ao)|C2eo will be important in the computation 0/5*1(0) in section\^ 

The previous result is a microlocal result and the subject of the present paper 
is to get a global result. 

In the adiabatic case (see Section EH]) , we get 



where gap(/x) is the minimal gap of the avoided crossing and A'_|_ are the slopes of 
the unperturbed eigenvalues at the crossing point. 

2 The local scattering matrix for the Landau- 
Zener normal form 

The goal of this section is to compute in a very explicit way the local 2x2 
scattering matrix T for the Landau-Zener normal form : 



70 (/i) 



gap(/i)2 



+ 0(/.^) 





(5) 
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with D = ^1-. 

I ax 

Let us put 7 = I a p and let us choose some small a > and assume 7 < . Let 
C^^ = {x^ = 7} be the characteristic manifold. The set C^^n{max(|x|, |^|) > a} 
is the union of 4 connected arcs. These arcs are labelled A™'°"* as follows: 

• A^^ = e I X < -a} 

. A-t = {(x,OGCLZ \ ^<-a} 
. A- = {(x,OeCLZ |^>a} 

• A°_^' = e I a; > a}. 

The meaning of the labels is as follows: 

• "in" (resp. "out") means that the arc oriented according to remark [2] is 
incoming (resp. outgoing). 

• "+" (resp. "— ") means that the vanishing eigenvalue of the dispersion 
matrix is the largest (resp. smallest) one. 

We start defining 4 WKB (exact) solutions of the previous system associated 
to the 4 Lagrangian arcs A™'°"*: 




Figure 2: the 4 arcs of the characteristic manifold and Morse indices 



X. 



x_ 



u^^'Hx) 



u'^{x) 



iX — l 



-1 



x) = ax_ 

5(0 = 
^?(0 = 



(6) 
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where /(^) is the /i— Fourier transform of f{x) defined by 



1 



f{0= / e-^^f{x)\dx\ , 

and x± = Y{±x)\x\ with Y the Heaviside function. 

Computing the Fourier transforms of u± and v±, we get the following com- 
patibility conditions in order to get microlocal solutions of ([5]) near the origin: 



If := y+iy;'^ + y.W'J" and PF""* := + are WKB-solutions 

of Equation (jS]) outside the origin, we get, for any microlocal solution near the 
origin. 



^+ J \ y- 

where T is the unitary matrix defined by: 
with 

i_^.^r(i + a) ^ i^.^r(i + a) 

A = ih-2+'l ^ Z_ ^^ e"^, B = ih-2^'i ^ Z_ ^^ e~"^ . (8) 
\/27ra v27ra 



The matrix T will be called the local scattering matrix associated to the 
singular point (and the choice of a normal form). Unitarity of T is checked using 
the well known formula 



nx 



sinh -Rx 



which implies \A\^ = \B\'^ + 1, AB = AB. 
The transmission coefficient 



1 

exp(-7r-) (9) 



gives the Landau-Zener formula. The previous explicit expression for the scat- 
tering matrix will allow to define in the next section the Bohr-Sommerfeld phases 
and to take into account interferences patterns due to several (avoided) crossings. 
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3 Singular Bohr-Sommerfeld phases 



3.1 Outline 

To each cycle c of the dispersion relation Co, we can associate, using the recipe 
of [9], a singular phase of the form 

S/i(/i) = So{fi) + hSi{fi) H 

where the S'/s are smooth w.r. to fi. 

In this section, we will define precisely these phases. We show that they are 
uniquely defined as formal power series w.r. to (/i, h). More precisely, each Sj{ii) 
is well defined modulo fiat functions on Z (see Remark [1]). We will give more 
precise properties of 5*0 in section H] and Si in section [5l So{fi) formal power 

series, a purely classical object derived from the dispersion relation, while 5*1(0) 
is a semi-classical object associated to phases given by the transport equation 
which in the adiabatic case are Berry phases. 

3.2 Bohr-Sommerfeld phases: a definition 

Let us take a simple oriented cycle c of the dispersion relation Co (boundary of 
a bounded component of T*R \ Cq). Let Zi, Z2, ■ ■ ■ , Zn he the singular points of 
c ordered cyclically around c. For each singular point Zj, let us build a FIO Uj 
and a "^DO gauge transform Aj (all depending smoothly on /i) which give the 
normal form of Theorem [1] with aj = aj (/i, h) a full symbol. 

We will define 'H^(c) = exp{iSh{l^) / h) as follows: we will denote by 
the images of W™' by the operators UjAj. These functions are WKB solutions 
of equation (jlj) associated to arcs of C^ near Zj. We introduce also WKB solutions 
Uj of (jl]) along arcs of C^ close to ]zj, Zj+i[. From those objects we get a global 
holonomy 'Hfj,{c) of the cycle c defined as follows: we have (by uniqueness, modulo 
multiplication by a full symbol, of WKB solutions) for example near Zj-. 

We define 'H^(c) = H"^ iVjX: . In other words, 'H^{c) is the holonomy of a sheaf 
on c given by the WKB solutions on the smooth part of the cycle and whose 
jumps of section are given from the normal forms. In our previous example 

Lemma 1 We have in^ic)] = 1 + 0{h°°). 
Proof. - 
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Figure 3: defining 7i^{c) 



Following [7j section 11.2.1. and Figure HI we associate to the cycle 
c an unitary scattering matrix Sc which is computable from the lo- 
cal unitary scattering matrices associated to the singular points and 
the holonomy 'H^(c). If this holonomy does not satisfy 17^^(0)1 = 
1 + 0{h°°), the global scattering matrix would not be unitary: the 
previous matrix is the product of (unitary) local scattering matrices 
and a diagonal matrix whose unique nonzero entry is l-L^^c). 




Figure 4: the scattering matrix associated to a cycle 



□ 

Taking the Logarithms, we get the phase Sh{lj)/h = Yl'jLo which is 

well defined modulo a multiple of 2it. 

Lemma 2 Given the gauge transforms used in the normal form of Theorem Ui 
the Taylor expansions of the Sj 's at fi = are uniquely defined. 

Proof. - 
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We will use the fact that the local scattering matrix computed in 
Section [2] is irreducible in the domain < < Vh, meaning 

that none of the entries are 0{h°°) in this domain. It implies that, 
for each j, the W^!^ °"*'s can be, up to a global multiplication by a 
symbol, defined as sets of WKB solutions for which the scattering 
matrix is given by Equations ([7]) and ([8]) with the value of a given by 
the normal form at the point Zj. 

□ 

4 The classical part Sq 

4.1 5*0 is classical 

We have the following: 

Theorem 2 For any simple cycle c of Cq, the function So{fi) depends only on 
the dispersion relation C^j,. 

Proof. - 

From the definition, 5*0 depends only on the terms in 1/h in the 
phases of the images by our normal form transformations of the ex- 
plicit solutions of the normal form. Those terms depends only on 
the canonical transformations used in the normal form and the as- 
sociated generating functions via stationnary phases (the Lagrangian 
manifolds) . 

□ 

4.2 5*0 as a regularized action integral 

Let us denote by = d{fi,Z). As in [9], it would be nice to get So{fi) as a 
regular isat ion of an usual action integral. A basic fact in [9j was that any simple 
cycle c is a limit of a cycle of as fi ^ 0±. This is no longer the case here 
because 7o > 0; one can see an example in section 1441 The idea is now to forget 
the initial problem and to work only with the dispersion relation which can 
be embedded into a larger family Ct^^ for which we can define action integral in 
some suitable sectors of the (t, /i) space. We can then restrict to t = and get 
our actions Sq. 

We will calculate Sq by first computing the same object So(t,/i) for Cj^^ = 
{p{x, ^, /i) — t = 0} and taking S'o(/i) = So(0, At). 

The cycle c is a limit of a cycle c(t, /x) of Ct^^ as {t, yu) — i- in some sector 
fi± := {(t,/i)| ±t > 0, « . 
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We have, for (t, /x) G fl, 



So 



(t,/i)= / + ^±7o,j(t,/i)(ln|7o,j(t,/i)| - 1) , 



where the ± signs depends on orientation and can be determined from the Loga- 
rithmic singularities of the action integrals. The contributions ±7oj(t, /i)(ln |7o,j(t, — 
1) come from the phase shift between x^^^^ and expressed as a WKB func- 

tion of the single variable x. 

Knowing that Sq is smooth, the previous formula defines the Taylor expansion 
of So w.r. to (t, /i) and hence the Taylor expansion of Sq w.r. to /i. 

4.3 The analytic case 

In the analytic case, we could also consider the Riemann surfaces = {p^ = 0} 
and look at some complex cycles on whose limit is c. Those cycles are not 
unique, but the real part of their action integrals are well defined and we can 
then take directly the previous regularisation. 

4.4 An example 

Let us consider the adiabatic equation: 



and the only cycle Cq of Cq passing by the singular points (±1,±1). It is clear 
that Co is not a limit of real cycle of C^, because the matrix has real 

eigenvalues for each t and so is the union of 2 disjoint graphs and has no real 
cycle. 

5 The subprincipal action 

We know that the Landau-Zener coefficient given by Equation is 0(/i°^) if 
>> \/h. It implies that in order to solve our problem up to 0{y/h) terms it 
is enough to know 5*0 mod 0(|/i|'^) and 5*1 for fi E Z. Let us assume that we have 
local coordinates so that E Z. We will describe below the calculus of 5*1(0). 

Lemma 3 Assuming \fi\ = 0, the principal part a{x)exp{iS{x)/h) of the WKB 
solutions of Equation ^ associated to arcs ]zj,Zj^i[ of Cq can be smoothly ex- 
tended beyond the singular vertices as WKB functions. 



hdX 
i dt 



A^{t)X 



with 
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Proof. - 



The property is invariant by FIO and it is enough to prove it for the 
solutions of the normal form given in Equation (j6]). The point is that 
ao(0) = 0, hence 7 = 0{h^). 

□ 

The previous result is related to the fact that the adiabatic theorem is still 
valid in case of eigenvalue crossings (see [HE]). 

We will define on c a piecewise smooth Hermitian line bundle L with a con- 
nection as follows: 

• On each arc [zjjZj+i], = ker ifciass(-2) with the connection given by the 
transport equation as in [12] (in the case of the adiabatic limit, it is the so 
called geometric connection or Berry phase [2]). 

• At each singular point, there are 2 limit fibers L±j and from Aq (defined 
in Remark [3|) we have an isomorphism between both limits given by trans- 
porting the isomorphism (1, 0) — )■ (0, 1) of C © on © C by (Ao)|c2eo- 

Definition 1 The phase exp(i5'7(0)) is the holonomy of the discontinuous line 
bundle L. 




m(c) = 2 - 3 X i 



Figure 5: the singular Maslov indices 

Using the calculus of [TU] (page 20) (we alert the reader that the previous 
convention for Maslov indices are not the same in the paper 0), we can also put 
the: 

Definition 2 The (singular) Maslov index m{c) G Z/2 of a simple cycle c which 
is the boundary of a bounded connected component o/T*M\Co is given by: m{c) = 
msmooth(c) + arising (c) whcrc msmooth(c) is thc usual Maslov index of a smooth 
deformation of c while msing(c) is a sum 0/ ±| associated to the singular points 
according to the rules of Figure O 
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The Maslov index of any cycle is defined by linearity from the previous Maslov 
indices, so it gives a cocycle. For example, the Maslov index of a smooth cycle 
(even if not simple) is the usual one, namely ±2. 



Theorem 3 Using the previous definitions, we have: 
Proof. - 

The proof follows essentially the lines of [2] P- 474-476. 

Let us give some details. A priori, there are several cases to check 
depending on the position of the cycle c at the singular points w.r. to 
the verticals. We will assume that the matrix 



x'io) = 




of the canonical transformation x = Xo satisfies a ^ and b 0, this 
is the generic case. We define 

_ f +1 if a6 > 
^ ~ I -1 if a6 < ■ 

The generating function (^(x, = ip2{x,y)+0{\x\'^+\y\'^) of x, defined 
X{y,-dy^p) = {x,d^ip), satisfies (a;,!/) = ^{dx^ -2xy + ay'^). We 
need to compute mod Oh(X) the values for x close to of the images 
by the normal form transform of 

and 

i^'M = ( ) . 

Let us assume that the principal symbol of the "^DO gauge trans- 
form is the N X N matrix 

(^{A){y,v) = { ^liy^v) a2{y,ri) ■■■). 

We get for the components of both WKB solutions for x small but 
nonzero: 

= (27r/i)-3/' / e^^^^'''y^+^y-y'^''^Cix,y)aiiy,ri)dydy'dri , 

Jy'>0 
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with C(0,0) = and 



Jri>0 

If we evaluate the integrals by stationnary phase, the dominant 
contributions come from the critical points and not from the bound- 
ary. The determinant of both Hessians are the same, while the signa- 
ture differs by 1. The final result follows then by 

• Looking at the value of the stationnary phase calculations as x 
is close to 0: the limits are respectively 

C(0,0)e^""/Vi(0,0) 

and 

C(0,0)a2(0,0) 

• if e > 0, one should add a contribution of the smoothed c, while 
if e < there is no such contribution. 

• Remembering that 

^o|c2eo = ( ai(0,0) a2(0,0) ) . 

□ 



6 Application 1: adiabatic limit with avoided 
crossings 

6.1 Adiabatic limit 



A,{hT)X (10) 



We consider the following equation: 

IdX 
i dr 

where A^{t), < t < a, is a self-adjoint matrix which is smooth w.r. to 

and we consider < r < a/h. The limit /i — > of this equation is called the 

adiabatic limit. 

We can rewrite Equation f fTOl) in a standard semi- classical form by puting 
t = hr: 

h dX . , s -.^ ^ N 

-^-Mt)x (11) 

where < t < a. 

We will assume that the eigenvalues of ^^(0) and A^(a) are all non degenerate. 
The scattering matrix S{n, h) : — t- is defined by X{0) — > X{a) where X 
is a solution of Equation ([11] ) . 
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6.2 Outside eigenvalues crossings 



Let X{t) be an eigenvalue of multiplicity 1 of Ao{t) for t in some open intervall /. 
Then Equation (JTT]1 admits a unique (up to multiplication by some function of 
h) formal WKB solution given by 



where A'(t) = \{t) and ao(t) satifies: 

• ao(t) e ker (Ao(t) - A(t)) 

• Vao(^) = where V is the geometric or Berry connection. 

Let us recall that V Q/Qta{t) = Uta'it) where Ut is the orthogonal projection 
of onto the eigenspace ker(Ao(t) — A(t)). 

The previous statement is the content of the so called quantum adiabatic 
theorem and goes back to [3j. 

6.3 Avoided crossings 

What happens when eigenvalues become degenerate at some values of t? 

Let us try to understand the generic situation. It is well known that eigenvalue 
crossings for a real symmetric (resp. complex Hermitian) matrix is a codimension 
2 (resp. 3) property. It is the content of the well known Wigner-Von Neumann 
theorem |21j . 

Physically, eigenvalues crossings can still occur for symmetry reasons. But, 
if we break the symmetry by a small perturbation of size n, we will get the so- 
called avoided crossings. We have now two small parameters: the semi-classical 
(adiabatic) parameter h and the perturbation parameter /i. The previous results 
allow to discuss the uniform expansion of the scattering matrix w.r. to both small 
parameters. 

6.4 Precise assumptions 

We will assume that the eigenvalues of Ao(t) cross transversally only by pairs on 
]0,a[. The dispersion relation C T*[0,a] is defined by p^{t,T) = det(rld — 
A^(t)). So that is exactly the union of the graphs of the eigenvalues of A^(t). 

6.5 Calculation of the scattering matrix 

Let us describe how to compute the global scattering matrix in the case of Figure 
[71 Let us start with the 4 local scattering matrices Sj, j = 1, ■ ■ ■ ,4 and the 2 
holonomies 'H^(cfc), A; = 1,2. 
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Figure 6: the dispersion relation for the adiabatic hmit 




Figure 7: recipe for the global scattering matrix 
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We try to describe a global solution of our system which is given from WKB 
solutions associated to each arc of a maximal tree of Cq. We have 10 equations 
with 13 unknowns which allow to compute y from x. 



w_ = H^{ci)w+ 




It turns out that the global scattering matrix is the product of 5 unitary matrices 
as follows : 




7 Application 2: EBK quantization rules 




Figure 8: the dispersion relation for the Born-Oppenheimer Hamiltonian 
We consider a Born-Oppenheimer Hamiltonian of the following form: 

where : M — > Herm(C^) is smooth w.r. to (x, v). We assume: 

17 



The eigenvalues of Vo(x) are of multiplicities at most 2 and cross transver- 
sally. 



• The following properness condition: 

V^{x) > p{x)ld 

where \im\x\_^oo p{x) = +00. 

• We choose E so that, for any x G M, -E is not a degenerate eigenvalue of 
Voix). 

• If the eigenvalue Xj{x) of Vo(x) satisfies Xj{xo) = E, then A^ (xo) 7^ 0. 

We can apply the previous method in order to compute EBK quantization 
rules for the equation {K^, — E)U = 0{h°°) . 

EBK quantization can be solved following the same path; but is this case we 
have the same number of equations than of unknowns and EBK rule is given by 
the vanishing of a suitable determinant as in |9] . 
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